Hierarchy of Higher Dimensional Integrable 
System 



YU Song-Ju t, Kouichi TODA t and Takeshi FUKUYAMA || 

Department of Physics, Ritsumeikan University, Kusatsu, Shiga 525 JAPAN 
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test. Nevertheless, it may be concluded that the BS hierarchy is integrable. 
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1. Introduction 

The inverse scattering transform(IST) method is a powerful tool for the investigation 
of (1 + 1) dimensional nonlinear differential equations. The extent of applicability 
of this method, however, depends on the dimension of space. The well known 
(2 + 1) dimensional integrable equations are the Kadomtsev-Pettviashivili(KP), 
Davey-Stewarson(DS) equation and some new integrable equations which has been 
constructed by B.G.Konopelchenko and V.G.Dubrovskyjl]. In this paper we consider 
a new integrable equation in (2+1) dimensions, the Bogoyavlenskii-Schiff equation (BS 
equation) |^ . This equation is given by 



where 



ut + ^{u)u^ = 0, (1) 



^u)^^dl+u+^u,d-' (2) 



is the recursion operator of KdV equation^. Potential form of this equation takes 
the form, 

0xt + j'Pxxxz + 4'x4>xz + ^4'xx4'z =0 (m EE (j)^). (3) 

The BS equation has an infinite number of conservation lawsj^, the Painleve 
property^, |^ and N soliton solutions^]. That is, the BS equation is the (2 + 1) 
dimensional integrable equation. For (1 + 1) dimensional case(2 = x), i.e. KdV 
equation, we can construct the higher order equations of the KdV(the KdV hierarchy) 
with Lax operators and the conserved quantities[^ ^. In this paper, we consider 
whether this construction of hierarchy is also valid for (2 + 1) dimensional integrable 
equation or not. Namely the purpose of the present paper is to construct the BS 
hierarchy and to check its integrability through the Painleve test. 

In Sec. 2, we review the hierarchy of an integrable equation in (1 + 1) dimensions 
and derive the Lax pair of the BS hierarchy. In Sec. 3, we investigate the integrability 
of the BS hierarchy obtained in the previous section. Sec. 4 is devoted to discussion 
mainly concerning with integrable equations in (3 + 1) dimensions. 

2. The Lax Representation of The BS Hierarchy 

The potential KdV equation 
1 3 

(I^Xt + ^042! + -^4'x(t>XX = 0, (4) 

admits the Lax representationpC[| 

[L,T]^0. (5) 
L, T operators have the form 

L{x,t) ^dl + Mx,t), (6) 
T{x,t)= (^L{x,t)i^ +dt=d,L{x,t)+T'{x,t)+dt, (7) 



where ( means the part of ( ) with non- negative powers of and 
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^Ax.-\4>^^d-^ + ---, (8) 

T'{x,t) = i^.a, - i^,,. (9) 

The second equahty in equation (7) is given for the later use of the BS equation and 
others. 

An infinite set of equations, 

[L,r„] = (j^=l,2,3,---), (10) 

constitute the potential KdV hierarchy and the equation for n is called the (2n+ l)th 
order potential KdV equation. The operator T„ has the form 

r„(x,<) = +9t =a,L(x,t)" + r:(a:,i)+9t. (11) 

For n ^ 2 and n — equation (|lO|) represents the 5th and the 7th order potential 
KdV equations, respectively: 



1 5/ . , 1 2 



l^xt -^ -^VQx -t- o I V^a; T^Vxx "t" VxVxx 



and 



(12) 



n - \^zdl - \^xxdl + Q(/.3. + \<i>l^ d., - ^(/.4. - ^C^xCpxx (13) 

for n — 2. 

<l>xt + '^'f'Sx + ^ ^50;^ + 10(j)x{(j)lx + (t>x4>Zx) + 202;05x + ^4>xx4>ix + 3032,^ = (14) 

and 



T' 

3 



\<l^xdl - ^(j>xxdt 



11 



11 

32 ' 



>2 . 

XX 



16^ ^ 19^3 ,o 

^<Pa^'?i3a; + Y^(Px ]Ox 



64 



32 



1 9 \ 

7 9 

!!'£!;'?!'4a: " -^(t>xx(t)?.x + ^( 



32 



>l(t>xx 



(15) 



for n = 3. Equations (12) and (|lj) are integrable[[ll|. 

The Bogoyavlenskii's method is to modify the form of T operator(7) in (1 + 1) 
dimensions to in (2 + 1) dimensions as 



r(a;, z, t) = d:,L{x, z, t) + dt + T'{x, z, t), 



(16) 
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where 

r{x,z,t) = ^(bA-\(i^..- (17) 

Namely we replace dxL in equation(7) by dzL. Here L has the same form as 
equation(6), though 0, in this case, is a function of x, z and t. Then, the Lax 
equation(||) becomes the BS equation We can construct the (2n + l)th order 

BS hierarchy analogously to the potential KdV hierarchy: 

[i,r„] = o, (18) 

where T„ have the forms 

T„(x,z,t) = 9,L(x,z,i)" + T,;(x,z,t) + 9t. (19) 

We could construct the same equations for the Lax equation(|l^ with other forms of 
the operator r„ (a;, z, i) = L ■ — L [dzL) L ■ — L +T^{x, z,t) + dt- Here the 

m — 1 times n — m times 
operators Tl^{x,z,t) is different from T^(x, z,t), though we have the same equation 
after all. For n — 2 and n = 3, equations(^8|) represents the 5th and the 7th order BS 
equations respectively: 
For n = 2, equation(n8h becomes 



4'xt ^" ^q4^xxxxxz ~^ i2^x4'xxxz ~t~ ^4^xx4^xxz i2^xxx4'xz ~t~ ^^xxxx^z ^x^xz 
+ l<l>x(l)xx<l>z + ^(l)xxdx\(f>l)z - 0, (20) 

and 



T2 = l<l>zdl - \cj)xzdl + {^Cjjxxz + \<l^x(l)z + \d~H(fi)^dx - ^(l)xxxz 

1 3 

- -j^iixiixz + ^(t)xx(t>z- (21) 

In this case, the above equation(20) is equivalent to that obtained by Bogoyavlenskiifi 



|12|. For n = 3, equation(18) takes the form 



H^xt ~t~ A T^xxxxxxxz ^ 'PxH^xxxxxz ~t~ t^^4^xxH^xxxxz ~^ ci'rxxx'Pxxxz 

64 16 32 8 

t~,rtr^XXXx4^XXZ ^ 4^XXXXx4^XZ ~^ r^XXXXXx4^ Z ~t~ , 4^x4^ Xx4^ XXZ 

61 lb 32 4 

+ '^'fixx'Pxz + ^4l4>xxxz + ^(Px^xxx^xz + Y^'t'x^xxxx^z + ^i'xx^xxxi'z 

11 15 

- -^4xxdx^{4>lx)z + -^(f'xxxxdx^i^Dz + 4>l4xz + -^(fi^xx^z 

+ l-4^^d-\(g)z + ^4x(l>xxd~\4x)z = 0, (22) 
lb lb 

and 
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1 3 \ 2 f ^ 5 

3 1 25 15 1 

+ g0a;0x£cz " ^^x^ ('/'^^a;). + Y^0a:a;a;0z + J^'l^l^^ + ^Q^x^^'t'^^^ 

1 2 \ ^ 3 1 3 

~t~ T7r0a:^x i^x^z ] (^x "^TT^xxxxxz ~r^4^x4^xxxz ~r^4^xx4^xxz ~^4^xxx4^xz 

16 / 64 16 16 8 

15 3 15 3 

+ ■:;7:0a;a:x2;0z - 70^0X2 + --;<j)x4'xx4>z + 7;7:4>xxd^^ {4>l.) z ■ (23) 

32 4 lb 32 

Equations(20) and (22) are, of course, reduced to the 5th and 7th order potential 
KdV equations in the case of z = a;, respectively. The calculations for n > 4 can be 
performed similarly, but the procedure to determine operators for larger n becomes 
more complicated. 

3. Application of The Painleve Test 

We have constructed the BS hierarchy in the previous section. In this section we 
proceed to check the integrability of the BS hierarchy (|l^) through the Painleve test 



in the sense of WTC|13, 14|| . For the BS hierarchy, there appears an ambiguity in this 
test. Hence we give the detail of calculation, which may help to understand what it is. 
In order to eliminate the integral operator, we rewrite the 5th order BS equation(20) 
in the form of coupled systems, 

113 1 

4^xt ^ ~r^4^xxxxxz ~^ 'Z^x4^xxxz ~t~ ~7 4^xx^xxz ~t~ 'z4^xxx^xz 

Id 2 4 2 

1 3 1 

+ ■^'t'xxxxi'z - Pxx4>xz + -^4>x<t'xx4'z - -^4>xxPxz = 0, (24) 

Pxx + 0' = 0. (25) 

We now consider a local Laurent expansion in the neighborhood of a non-characteristic 
singular manifold 7(2;, z,i) = 0, (7. 7^ 0). Let us assume that the leading orders of 
the solutions of cquations(24) and (25) have the forms 

= 007", P = Pol''- (26) 



Here 0o and po are some analytic functions. Substituting (|2^) into (24) and (25), and 
equating the powers of the most dominant terms, we obtain 

a = -l, /3 = -2 (27) 

with 

0o(0o-67,)(0o-27,) = O, 6po + 0o-O. (28) 
To find the resonances we now substitute the full Laurent expansion of the solutions 
= ^0^y-i, p = ^p^Y-2, (29) 

j=0 J=0 
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into equations(24) and (25). Rearranging equation(24) into terms of 7^ ^ and the 
other higher powers of 7 (fj), we obtain recurrence relations for (f>j, pj, 



[U - m - 2)(j - 3)(j - 4)(j - 5)(j - 6)7^7. 

- W{j - - 6){f - 7 J + 16)7^7.00 

- 108(j - - 2)737,^0 + 12(j - l)(j - 6hh.<Pl] h 



36(j - 2)(j - 3)7^7.0oP,)7^"' = f, 



(30) 



Similarly rearranging equation(25) into terms of 7^ ^ and higher powers of 7 {gj), we 
have 

(-2(i - + (j - 2)(j- - 3)7|p,)7^-^ = .9,. (31) 

Here Jj and are given in terms of and p; (0 < Z < j — 1). Then we get three 
types of resonances: 



j = 0,1, 1,2, 3, 4, 5, 6, (32) 
J = -1,1,2,2,3,5,6,8 (33) 

and 

J = -3,-1,1,2,3,6,8,10. (34) 
These correspond to the choices of solutions of equations (|28 



and 



= 0, po - 0, (35) 
60 = 27,, po = -(2/3)7^ (36) 



67^, Po = -67x, (37) 



respectively. Here the first choice (35) does not include the universal resonance j = —1 
and is discarded. If the recurrence relations are consistently satisfied at the resonances 
then the differential equations is said to possess the Painleve property. 

Succeeding coefficients and pj are determined from equations (30) and (|3l|). 
However, from the consistency condition, they must include arbitrary functions at 
the resonances. To simplify the calculations, we use the reduced manifold ansatz of 



Kruskalj 15 . That is, 7(2;, z,t) = x + 5{z^t) and 0j , pj are function of z and t. 

Equations (30) and ( ^ ) must be satisfied in the respective powers of 7. The lowest 
powers of 7, (7~'',7~'*) in equations ((30), (|3l|)) give 

0o(0o-2)(0o-6) = O, po = -^0o, (38) 
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which constitutes a set of equations(35), (36), and ( pT\). Higher powers of 7,(7 

with positive integer k) in equations ( ( 30 ) , (pi]) ) lead us to the foUowing 



7" 



-4+k 



consistency conditions. 

(7"', 7"'): (boPiS, = 0, pi=0. 

Hence 0i is arbitrary. 



(7-^7-^) 



tiQ ~ 2)0i,z + 2(00 - 5)(/)2(5z 



= 0, 



= 0, 



(39) 



(40) 



These equations imply that 0o is fixed as 2 and p2 is arbitrary, which contradicts 
equation(p7|). So we must choose (33). 

(41) 
(42) 



(7"', 7"') : (300 - 20)03^. = 0, 03 - 0, 
and p3 is arbitrary. 

(7-3, 7O) : 6(1304 - 3p4)S, + 86t - P3,z = 0, p4 = 604- 



(7-^ 7^) : 9(805 - SP5)S. + (604 - P4). = 0, ps = 7:^5 



and 



(43) 



(7-1,7'): 36(506 - 3p6)5. + (805 - 3p5). = 0, P6=306. 



(44) 



It follows from equations (^) and (^4|) that one of the two variables must be arbitrary 
in both pairs (05, p^) and (06, pe). 



(7°, 7') 



12(1307 - I5p7)d, - (506 + 4p6). 



hcpl.z = 0, P7 = ^<A7- (45) 
5 



(7\7') 



(50408 - 540p8) + 720^ 



+ 304(-18p4'^z + SSt - P3,z) 



P8 



+ 2(607 



5p7)z = 0, 

10^^ 



(46) 



That is, one of the two variables (08, ps) must be arbitrary. Problematic is the 
following point: 02 at double root j — 2 has been fixed as and the system lacks one 
more arbitrary function. However, 0^- and pj are determined consistently and, thus, 
dependent variables and p have only movable poles. The apparently similar situation 
occurred in the system of coupled non-liner Schodinger equation(NLS) ]l6t . Starting 
from NLS leaving the constant coefficients as free parameters in ||l^, Sahadevan et 
al. these parameters so as to include N arbitrary functions at iV-tuply degenerate 
resonance. However, there seems to exist an essential difference between their NLS and 
our 5th order BS equation since NLS does not allow the consistent Laurent expansion 
except for the special choices of constant coefficients above mentioned. In the 5th order 
BS equation there is no such room to adjust the number of arbitrary function, though 
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the 5th order BS equation has no inconsistency. We want to add one remark here. 
The coupled system of equations (24) and (25) is reduced to the 5th order potential 
KdV equation in the case of z = x. Even in that case, the resonance structures (the 
positions and the number of arbitrary functions etc.) are quite the same as the 5th 
order BS equation. The 5th order potential KdV equation is manifestly integrable. 
From these facts we may conclude that the 5th order BS equation is also integrable. 

Analogously, we have checked that the 7th order BS equation is integrable in the 
same sense. 

Higher order equation can be constructed in several different ways. All the higher 
order BS equations(|l^) are rewritten as 

ut + $(«)"«, = (47) 

by the use of the recursion operator of KdV equation(^). Equation(|47|) was conjectured 
to be integrable because the integrability arose from the existence of the recursion 
operator[|l^. However, the recursion operator in (2 + 1) dimensions is not so clear as 
in the case of (1 + 1) dimensions. Explicit calculation in this section has uncovered a 
new feature and has shown that the BS hierarchy is integrable. 



4. Discussion 



In this paper, we have constructed the BS hierarchy (p^ using the Bogoyavlenskii's 
method. This hierarchy is represented by the recursion operator (^) and reduced to 
the potential KdV hierarchy by setting z — x. The BS hierarchy(18[), ( J^f\ j has been 
proved to be integrable through the Painleve test. The BS equation is an extension of 
the KdV equation like the KP equation is. However their extensions differ in algebraic 
structure, which is schematically depicted in Figure.|l|. 



KdV 



(48) 



di) 



KP 



BS 



XI 



X3 



Figure 1. Scheme of extensions of the KdV equation. There are two directional 
routes of extensions: One leads us to the BS equation and the other does to the KP 
equation. X equations are possible (3 + 1) dimensional integrable equations obtained 
from further extensions. 

The potential KP equation is obtained from the potential KdV equation by the 
replacement. 



L{x,t) — > L'{x,y,t) =dl + Mx,y,t)+dy 
EE L{x,y,t) +dy, 

The form of T{x, y, t) is 

T(x, y, t) = a,L(x, y, t) + T\x, y, t) + dt, 



(48) 
(49) 
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where 

1 13 

T'{X, y, t) = -(j),;d.j; - -(j),;^ " - ^j, . (50) 

Here we have denoted an extended spatial coordinate as y. These extensions are 
contrasted with those of the BS equation(|l^) in which explicit y (in that case z) 
dependence has appeared not in L but in T. 

These two extensions may be performed in the different spatial dimensions (i.e. 
y z), which is described by XI in Figure.^ For instance, we consider the Lax 
equation(||) with the operators L' and T having the form 



where 



L'{x,y,z,t) ^dl + (l)x{x,y,z,t) + dy ^ L{x,y,z,t) + dy, (51) 
T{x, y, z, t) = dzL{x, y, z, t) + T'{x, y, z, t) + dt, (52) 



T'{x, y, z, t) = - -6.,, - -9- Vy. + cdl (53) 



and c is a constant. However, the Lax equation(||) with operators(51) and (52) in 
(3 + 1) dimensions is reduced to (2 + 1) dimensional equation, 

0xt + -^(Pxxxz + i'xi'xz + -j^i'xxi'z + C^d^^(f)zzz = (54) 

with the condition, 

z = 2cy. (55) 

This condition comes from the requirement that the coefficient of dz should be vanished 
in the Lax equation and, therefore, seems to be indispensable. Equation (^^ in (2 + 1) 
dimensions is the same as that obtained by BogoyavlenskiiQ. We have checked that 
equation(|5^) passes the Painleve test. By the dependent variable transformation 
4> — 2(logr)a;, equation ( p4[ ) is transformed into the trilinear form, 

(36Ti^7i + T,%* + 8T,'%*% + 9T/)r • r ■ r = 0. (56) 

The operators T, T* are defined by|||, |l9| 

%"f{z)-g{z) ■ h{z) = {dz, +jdz, +fdz,rf{z,)g{z2)h{z3)\z,=z,=zs=z, (57) 

where j is the cubic root of unity, j — exp(2i7r/3). T* is the complex conjugate 
operator of % obtained by replacing {dz, + jdz^ + fdz^) by {dz, + fdz^ + jdz^)- 
Equation(pq) coincides with the trilinear form coming from the quite different 
approach Ml. Namely, Hietarinta, Grammaticos and Ramani derived this equation 
from the singularity analysis of trilinear equations. We have obtained N soliton 
solutions of equation (pq) by the computer program Mathematica up to A'^ = 5. 



However, we can not construct the hierarchy of equation(54). Namely we can not 
determine T'{x, y, z, t) operator in 



T{x, y, z, t) = dzL{x, y, z, t)" + T'{x, y, z, t) + dt, 



(58) 
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so as to constitute Lax equation(|^). 

A little bit different formulation to reach XI candidate was proposed by us|^, 



Ut + ) + -Uyy = 0. (59) 



Unfortunately, equation (p9| ) has a movable logarithmic branch point. 

One of the examples for (3 + 1) dimensional X3 equations in Figure.^ is obtained 
also by Bogoyavlenskii|P. However, it is not affirmative unless it is proved to be 
integrable. 

The respective routes in Figure.|l| are not unique. Different formulations at each 
path are possible, though we can not specify which is correct at the present stage. 
Such specification enforces us to clarify the algebraic structures and integrability of 
the equations depicted in Figure. |l| and their hierarchies furthermore. 
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